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Abstract. - We present a new method to compute the electronic structure of correlated materials 
combining the hybrid functional method with the dynamical mean-field theory. As a test example 
of the method we study cerium sesquioxide, a strongly correlated Mott-band insulator. The 
hybrid functional part improves the magnitude of the pd-band gap which is underestimated in the 
standard approximations to density functional theory while the dynamical mean-field theory part 
splits the 4/-electron spectra into a lower and an upper Hubbard band. 



Introduction. — Recently, there has been consider- 
able progress in the realistic description of strongly corre- 
lated materials by combining density functional theory [1] 
(DFT) with the dynamical mean-field theory (DMFT) 
[2-5]. In this DFT+DMFT approach [6,7], DFT is em- 
ployed to obtain an effective mean-field description of the 
weakly correlated bands while the local correlations in the 
more strongly correlated bands (i.e. the d-bands of transi- 
tion metals and the /-bands of Lanthanidcs and Actinides) 
are treated exactly. The DFT+DMFT method allows to 
predict spectra and energies of stronlgly correlated mate- 
rials. 

The DFT+DMFT method has been successfully applied 
to a variety of interesting materials that conventional band 
structure theory is unable to deal with. For example, 
using the DFT+DMFT method the 25% of volume in- 
crease in the transition from the a- to the <5-phase of Pu 
could be explained by the presence of strong correlations 
in S-Pvl [8]. However, by construction the DFT+DMFT 
approach, does not work so well in situations where the 
one-electron spectra of the weakly correlated bands are 
not well approximated by the Kohn-Sham (KS) spec- 
tra of DFT. For example, the local density approxima- 
tion (LDA) and the generalized gradient approximation 
(GGA) notoriously underestimate the magnitudes of band 
gaps of insulating materials. On the other hand within 
the chemistry community, very accurate functionals called 
hybrid functionals [9], have been constructed by mixing 
LDA/GGA functionals with Hartree-Fock. The hybrid 
functional (HYF) approach has been tremendously suc- 



cessful in providing very accurate energies for molecules. 
Moreover, one-electron spectra computed with HYFs give 
fairly accurate gaps for semiconducting materials [10]. 

In this work we propose a new method that combines 
the HYF approach with DMFT (HYF+DMFT) to yield 
a quantitatively and qualitatively correct description of 
combined band and Mott-Hubbard insulators. The HYF 
part improves the effective static mean-field description of 
the uncorrelated electrons while the DMFT part describes 
the dynamical local electronic correlations of the strongly 
localized electrons that can lead to Mott insulating be- 
haviour. Although the HYF approach introduces a new 
parameter, a that determines the amount of Hartree-Fock 
exchange, we show that this a can actually be linked to 
the Coulomb repulsion parameter U of the DMFT calcu- 
lation. 

An important example which illustrates the need for 
the HYF+DMFT method is provided by the rare earth 
scsquioxides [11] series which are insulators. In addition 
to the Mott-Hubbard gap between the occupied and the 
unoccupied 4/-bands, a band gap between the uncorre- 
lated O 2p- and 5d-bands opens. One thus has to deal 
with the two-fold problem of finding an accurate descrip- 
tion for the pd-band gap and the 4/-Hubbard bands in the 
same material. 

DFT and related static mean-field methods fail to de- 
scribe the splitting of the 4/ Hubbard band, without in- 
voking some form of magnetic long range order. For ex- 
ample, HYFs give the correct magnitude for the band gap 
of antiferromagnetic (AF) Ce203 [12,13]. However, the 
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symmetry-breaking is essential to obtain the correct insu- 
lating behaviour although it is clear that the magnetism 
should not be the driving mechanism behind the insulat- 
ing behaviour since the Nccl temperature is only 9 K and 
thus much smaller than the measured gap of about 2.5 eV. 
Similarly, symmetry-breaking is crucial in order to capture 
the insulating behaviour of Ce2C>3 and of Ce2S3 with the 
DFT+U method [14-17]. Thus a satisfactory description 
of CC2O3 should also yield the correct insulating behaviour 
in the paramagnetic phase. 

Correlations on top of DFT, can be added within the 
DFT+DMFT method. Indeed recent studies of the rare 
earth oxides have been carried out within this approach. 
The DFT+DMFT [18] approach, succesfully describes the 
opening of the Hubbard band in the 4/-shell of these ma- 
terials. However, they underestimate the pd-band gap. It 
is natural then to combine the virtues of the HYF and the 
DMFT method, in a HYF+DMFT approach, which is the 
subject of this paper. 

Method. - The HYF+DMFT method follows the 
well established DFT+DMFT methodology. We focus on 
the one-particle Green's function: 

G(k,Lu)^(u; + ti-H(k)-t(uj))- 1 (1) 

which is expressed in terms of a one-body Hamiltonian 
H(k) and a local self-energy E(w). The self-energy de- 
scribes the dynamic electron correlations of the strongly 
localized electrons and thus has non-zero elements only in 
the block of correlated bands. The correlations are cap- 
tured by a Hubbard-like term which is added to the one- 
body Hamiltonian H(k) in the strongly-correlated sub- 
space: 

7~tu = 2 y ] U ai a 2 b 1 b 2 Ca 1(Tl cj l2(T2 Cb 2( j 2 Cb lCri (2) 

where the indices 01,02,61,02 denote orbitals of the cor- 
related subspace in a local basis set. 

Within DFT+DMFT the one-body Hamiltonian H(k) 
is given by the effective KS Hamiltonian Hk s — — <^;V 2 + 
^ext( r ) + Wh(r) + v xc (r). The KS Hamiltonian gives an 
effective one-body description of the electronic structure 
taking into account the Coulomb interaction on a mean- 
field level by the Hartree potential Wh and the exchange- 
correlation (XC) potential v xc - 

One of the main shortcomings of the standard approx- 
imations to DFT like LDA and GGA is their difficulty to 
describe insulating materials. On the one hand the insuf- 
ficient cancellation of the self-interaction error by the ap- 
proximate LDA and GGA XC functionals results in band 
gaps that are generally too small compared to the mea- 
sured band gaps in semiconductors and insulators [10]. 
On the other hand DFT is strictly speaking only a ground 
state theory. Thus it does not necessarily give a correct 
description of excited state properties like the band gap of 
insulating materials. 



The HYF method [9] improves on LDA/GGA by intro- 
ducing a fraction a of exact Hartree-Fock exchange into 
the XC potential: 

= av^ + (l-a)v^ + v^ (3) 

where w| ga is the LDA/GGA exchange potential and £f ga 
is the LDA/GGA correlation potential, v^f is the Hartree- 
Fock exchange potential which is a non-local (i.e. non- 
diagonal in real space) effective one-body potential: 

^ f (r,r') = <r| |r'> = -^(r, r')V ee (r - r') (4) 

where p(r, r') is the density matrix and V ee (r — r') = 
(r, r'| V ee r, r') = l/||r — r'|| is the bare Coulomb inter- 
action between two electrons. 

One can think loosely of the HYF approach as a first 
order correction of the KS Hamiltonian in a fraction a of 
the bare Coulomb interaction V ee . This aV ee can be inter- 
preted as a "screened Coulomb interaction" in a similar 
way as the [/-parameter in the LDA+U method [14, 19]. 
But in contrast to the LDA+U method, HYFs make the 
correction in the screenend Coulomb interaction on the en- 
tire LDA/GGA Hamiltonian and not only within a small 
subspace of atomic orbitals. Also note that the XC po- 
tential of the HYF method is non-local (i.e. dependent on 
both r and r' instead of r alone) due to the contribution 
of exact Hartree-Fock exchange to the HYF. This shows 
that the HYF approach really falls outside the framework 
of conventional DFT where the XC potential is required 
to be local. 

It turns out that the optimal amount a of Hartree-Fock 
exchange is almost universally of 20%-25% which has been 
rationalized by perturbation theory [20] . In the following 
we employ a functional similar to the popular B3LYP and 
B3PW functionals of quantum chemistry [21] that mix 
20% of Hartree-Fock exchange with 80% of GGA exchange 
functionals but instead of GGA we employ plain LDA here 
and refer to this HYF as LDA20. 

By adding the Hubbard term (2) to the KS band struc- 
ture some of the Coulomb interaction within the correlated 
subspace is double-counted since it has already been taken 
into account on a static mean-field level by the KS Hamil- 
tonian i/ks- Thus the KS Hamiltonian has to be corrected 
by a double-counting correction (DCC) term: 

H(k) = H ks (k) - h dc (5) 

In the case of LDA/GGA an exact expression for this 
DCC term h dc is not known, and several forms of the 
DCC have been suggested [22]. In the case of HYFs, how- 
ever, at least the DCC for the Hartree-Fock contribution 
to the HYF is known exactly: It is the Hartree potential 
■Oh plus the Hartree-Fock exchange potential projected 
onto the correlated subspace. Below we argue that this 
is the only relevant contribution to the DCC term of the 
HYF+DMFT approach. However, we cannot proof this 
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conjecture strictly so that ultimately it must be justified 
by its success in the application to different materials. 

Our argument goes as follows: Bearing in mind that the 
HYF approach corresponds to a first order correction of 
the effective LDA Hamiltonian in a fraction a of the bare 
Coulomb interaction V ee , one first has to remove entirely 
this perturbation correction in aV ee in the correlated sub- 
space since the DMFT calculation will treat the Coulomb 
interaction (locally) exact in that subspacc. Thus if we 
take aVabcd = a (ab\ V ee \cd) as the screened Coulomb in- 
teraction Uabcd in the correlated subspace we obtain for 
the DCC of the HYF approach: 

(a<j\ h*g \ba) = a (aa\ v h + \ba) + (aa\ ti^ c \ba) (6) 

where the last term (aa\ h}^ c \b, a) corrects the double- 
counting in the LDA correlation potential u' da within the 
correlated subspace. This term is not known exactly since 
the contribution of a subspace to the total LDA corre- 
lation potential cannot be calculated exactly. However, 
the correlation potential is usually much smaller than the 
exchange potential, and we will thus neglect this contri- 
bution here. 

Assuming that only the direct Coulomb interactions 
U = Uatab and the exchange Coulomb interactions J = 
Uabba arc important, we find the following simplified term 
for the on-site HYF DCC which allows us to predict the 
positions of the strongly correlated orbitals for the DMFT 
calculation: 

(aa\ \aa) « U(N f - <) - J(NJ - <) (7) 

where Nf is the total number of correlated electrons per 
atom and is the number of electrons in atomic orbital a 
with spin a. Eq. (7) is different from the usual expression 
for the LDA DCC in that the HYF DCC term becomes 
now orbital dependent so that unoccupied orbitals experi- 
ence a larger shift than occupied ones. Also note that since 
the Hartree-Fock potential is non-local (i.e. non-diagonal 
in real space) the DCC term is also non-local. But by 
construction the DCC term only acts on the correlated 
subspace. 

Within DMFT the self-energy is determined self- 
consistcntly by mapping the original problem onto an An- 
derson impurity problem. To this end the local Green's 
function, 

G loc H = ^ k (w + M -i?(k)-SH)- 1 (8) 

projected onto the correlated subspace is equated to the 
Green's function of the equivalent impurity problem: 

GoM = (^ + m-^o-AH-E / ( w ))- 1 (9) 

Here Ho is the (single-particle) Hamiltonian of the im- 
purity site, A(w) is the hybridization function with the 
conduction bath electrons, and S/ = P/EP/ is the full 
self-energy projected onto the correlated subspace where 



Pf is the projection operator for the correlated subspace. 
The mapping Gq(uj) = PfGi oc (oj)Pf defines the so-called 
self-consistency condition which is central to the DMFT 
method. The mapping yields the hybridzation function, 

A(w) = w + ft - H - t f (u) - (PfGio^Pf)- 1 (10) 

with Ho = Ylk PfHityPf. These are the relevant quanti- 
ties for solving the impurity problem. By solving the im- 
purity problem, one obtains in turn the self-energy S(w). 
Equations (8)-(10) define the self-consistent DMFT pro- 
cedure for computing the self-energy S(w). 

Solving the impurity problem given by (9) is the com- 
putationally most demanding step in most DMFT calcu- 
lations. A variety of solvers — each suitable for a certain 
region of parameters — is available to deal with the impu- 
rity problem. An overview over the different techniques 
can be found in e.g. Rcf. [5]. Since we are interested in 
describing a Mott insulator a suitable method for solving 
the impurity problem is an expansion in the hybridiza- 
tion strength in the so called non-crossing approximation 
(NCA) [23]. 

Results. — In order to show how the HYF method 
efficiently improves the gap of band insulators in compari- 
son with conventional DFT methods we first perform LDA 
and HYF calculations of La203 which is similar to Ce203 
but does not have the strongly correlated 4/-electrons. It 
is a typical band insulator with a reported band gap of 
about 5.5 eV between the O 2p-valence band and the La 
5c?-conduction band [11]. For the LDA and HYF calcula- 
tions we employ here and in the following LDA+DMFT 
and HYF+DMFT calculations of Ce 2 3 the CRYSTAL06 
ab initio electronic structure program for crystalline solids 
[24] together with a Gaussian basis set and pseudo poten- 
tial by Cundari and Stevens tailored for the rare earth 
series elements [25]. The basis set retains the 4/, 5s, 5p, 
5d, 6s and 6p shells for the electronic structure calcula- 
tions while the inner shells arc described by an effective 
core pseudo potential. In the self-consistent LDA and 
HYF calculations the irreducible Brillouin zone is sam- 
pled by 193 fc-points while in the DMFT calculations the 
entire Brillouin zone is sampled by 1000 /c-points. The 
DMFT calculations are performed at a finite temperature 
of k B T = 0.2 eV. The LDA+DMFT and HYF+DMFT 
calculations are not fully self-consistent in the charge den- 
sity [18]. 

Fig. 1 shows a comparison between the density of states 
(DOS) calculated with a pure LDA functional on the one 
hand and the LDA20 HYF on the other hand. While 
in the LDA calculation the magnitude of the band gap 
is underestimated as expected by more than 1.5 eV, the 
HYF calculation does indeed give the experimental value 
of the band gap of about 5.5 eV. 

We would like to stress here that changing the LDA or 
GGA functional part of the HYF does not alter the results 
significantly. In fact, we have obtained very similar re- 
sults and the correct magnitude of the band gap also with 
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Fig. 1: (Color online) DOS of La203 calculated with an LDA 
functional (dashed red line) and with the LDA20 HYF (con- 
tinuous blue line). 

the popular B3LYP or B3PW functionals. In contrast 
the results depend quite strongly on the exact amount of 
Hartree-Fock exchange. Thus increasing the amount of 
Hartree-Fock exchange to 25% already increases the mag- 
nitude of the band gap to about 6 eV. 

To demonstrate the above developed HYF+DMFT 
method we apply it now to the exemplary case of Ce 2 03. 
As before we employ the LDA20 HYF which reproduces 
correctly the band gap of La2C>3. The starting point is 
a HYF calculation of the paramagnetic phase of Ce2C>3 
which results in a metallic state with the Fermi level in 
the 4/-band. In order to prevent excessive symmetry 
breaking of the Ce 4/-orbitals prior to the DMFT cal- 
culations, the calculation is done at a finite temperature 
of k B T w 0.5 cV. 

In spite of the quite high temperature and although 
the crystal field splitting is actually quite weak (of or- 
der 0.1 eV), the orbital symmetry of the 4/-orbitals has 
been broken: The energy difference between lowest and 
highest 4/-orbital is about 1.2 cV. The reason for the 
symmetry-breaking is the Hartree-Fock contribution to 
the HYF which tends to break symmetries to lower the 
energy of the system so that those 4/-orbitals that are 
slightly favored by the weak crystal field splitting become 
more occupied during the self-consistent solution of the 
KS equations while those unfavored by the crystal field 
splitting become less occupied. Therefore the splitting 
is strongly enhanced by the Hartree-Fock term. Since 
the Hartree-Fock term is the principal responsible for the 
symmetry-breaking the energy difference between lowest 
and highest 4/-orbital is reduced to 0.4 cV when applying 
the HYF DCC scheme outlined above, partially restor- 
ing the symmetry This can be understood by consider- 
ing the orbital dependence of the HYF DCC (7) which 
shifts occupied orbitals less than unoccupied orbitals de- 
creasing the symmetry-breaking. At lower temperatures 
the symmetry-breaking becomes even stronger: E.g. for 
fcgT rs 0.2 eV the splitting of the 4/-orbitals is still of 
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Fig. 2: (Color online) Tot. DOS (dashed red line) and partial 
Ce 4/ DOS (continuous blue line) of Ce2 03 calculated with 
LDA+DMFT (a) compared to the HYF+DMFT approach (b) 
explained in the text. The on-site Coulomb repulsion for the 
DMFT calculation is U — 5 eV and J = 0.2 eV in both cases. 

about 2.3 eV after applying the DCC. 

We note that the on-site Coulomb repulsion aV aa aa cor- 
responding to the Hartree-Fock contribution is about 5 eV. 
This value is only slightly smaller than the U usually em- 
ployed for the Ce 4/-orbitals in actual LDA+U an DMFT 
calculations [15,16,18] which is between 5.5 eV and 6.5 eV. 
This again points to the correctness of the interpretation 
of aVaaaa as a screened Coulomb interaction similar to the 
U of the LDA+U method [14]. 

In Fig. 2 we compare DMFT calculations (a) on top of 
the plain LDA calculation and (b) on top of the LDA20 
HYF calculations. In both cases the 4/-band splits into 
a lower Hubbard band filled with one electron per Ce 
atom and an upper unoccupied Hubbard band. We take 
ctVaaaa = 5 eV as an estimate for U and similarly aV a bba = 
0.2 eV as an estimate for J. The resulting band gap be- 
tween the occupied 4/-Hubbard band and the empty 5d- 
conduction band is about 2.3 eV in both cases which is 
in very good agreement with the measured band gap of 
2.5 eV [11]. The agreement can be improved by slightly 
increasing U to about 5.5 eV. Most importantly, the pd- 
gap between the occupied O 2p-bands and the unoccupied 
Ce 5d-bands is 5.5 eV in the HYF+DMFT calculation 
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and is thus exactly the experimentally measured pd-gap 
of about 5.5 eV. This is a considerable improvement over 
the LDA+DMFT where the pd-g&p is only of about 4 eV. 
The occupancy Nf of the Ce 4/-shcll is very similar in 
both methods: We find Nf = 1.00 for the LDA+DMFT, 
and N f = 0.98 for the HYF+DMFT calculation. Finally, 
we mention that we have also obtained satisfactory re- 
sults (not shown) for the Nd2 03 compound using the same 
methodology 

Conclusions. — In conclusion, we have proposed 
a new method that combines the HYF approach with 
DMFT. We have shown that this HYF+DMFT method 
gives a qualitatively and quantitatively correct descrip- 
tion of the electronic structure of Ce2 03 as a prototyp- 
ical example of the strongly correlated insulating Rare 
Earth compounds. While the HYF part fixes the mag- 
nitude of the pd-band gap which is underestimated in 
the LDA+DMFT approach, the DMFT part takes care of 
the strongly localized 4/-electrons which are not properly 
taken into account within conventional KS band theory. It 
thus predicts correctly the opening of the Mott-Hubbard 
gap in the Ce 4/-band in addition to the band gap between 
the O 2p-valence band and the Ce 5d-conduction band. 
While the HYF is computationally slightly more expensive 
than the LDA, its cost is either comparable or less than 
that of the DMFT calculation. It is also computationally 
considerably less demanding than the GW approximation 
which gives an accurate description of the quasi-particle 
spectra of weakly correlated materials [26, 27] and can be 
combined with DMFT in a natural way [28]. Hence, the 
HYF+DMFT approach is an attractive avenue for improv- 
ing the overall accuracy of spectra in materials containing 
both correlated and uncorrelatcd electrons at reasonable 
computational cost. 
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